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Rare flavour changing neutral current decays are sensitive indirect probes for new effects 
beyond the Standard Model (SM). In the SM, these decays are forbidden at tree level and are 
therefore loop-suppressed. In SM extensions, new, heavy particles can significantly contribute 
and affect both their branching fractions as well as their angular distributions. 

The rare decay —>■ K'^'K~)pUpL~ is of particular interest, since it gives access to 

many angular observables, allowing to model-independentiy test the operator structure of the 
decay. A previous analysis of the angular distributions of the final state particles showed 
interesting tensions with SM predictions using the data sample taken by the LHCb detector 
during 2011. These proceedings will summarize latest results on rare decays from the LHCb 
experiment with emphasis on the angular analysis of the decay —>■ K*^pUpT, using the 
full Run I data sample of the LHCb experiment. 


1 Introduction 

Rare flavour changing neutral current (FCNC) decays are, in the Standard Model (SM), forbid¬ 
den at lowest perturbative order and proceed via loop-order diagrams. New heavy particles in 
extensions of the SM can appear in competing Feynman diagrams and significantly affect both 
the branching fractions of rare decays and the angular distributions of the final state particles. 
Studies of rare decays therefore constitute sensitive searches for effects beyond the SM, and 
furthermore allow to probe the underlying operator structure via global fits^’^’^’^. 

2 Angular analysis of the rare decay ^ 

2.1 Angular observables in )■ ii~ 

The rare decay —)■ is of particular interest since the A+vr p final state allows 

access to many angular observables. The final state is fully defined by the three decay angles 
= {cos 6 i,cos 6 K,(f>), and the invariant mass of the dilepton system squared^. The CP- 
averaged angular distribution of the decay B^ ^ in a bin of is given by 
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Figure 1 - (Left) Invariant mass of the K'^n~fi~ system vs. q^. The signal decay —>■ is clearly 

visible as vertical band. (Middle) The signal decay ^ integrated over q^. The signal yield is 

2 398 ± 57. (Right) The high-statistics control mode >■ 


Here, Fl denotes the longitudinal polarization fraction of the K*^ and ^fb the forward-backward 
asymmetry of the dimuon system. The LHCb collaboration performed two angular analyses 
using angular folding techniques to determine the observables with the data taken during 2011, 
which corresponds to an integrated luminosity of lfb~^. While all the angular observables in 
Ref. ® are found to be in good agreement with SM predictions, the measurement of the less 
form-factor dependent observable = S^/ \/Fl(1 — Fl), that was proposed in Ref. ®, shows a 
local deviation from the SM prediction, corresponding to 3.7 standard deviations (cr) Global 
hts show that the tension can be reduced by a negative shift of the Wilson coefficient Cg which 
parametrises the vector coupling strength New physics (NP) explanations for this shift 

include the possibility of heavy Z' gauge bosons®d0di.i2,i3 leptoquarks^^d5,i6,i7^ However, the 
signihcance of the tension could be reduced if hadronic uncertainties are underestimated 

2.2 Selection of signal candidates 

An update of the angular analysis of —)• using the full Run I data sample corre¬ 

sponding to 3fb~^ was eagerly awaited in the community and preliminary results are presented 
here for the first time^^. The selection of —?■ K*^gL~ signal candidates is improved compared 

to Ref. with a simplihed, yet more efficient, multivariate classifier to reduce combinatorial 
background events and more stringent vetoes to reject peaking backgrounds. Figure 1 gives the 
distribution of the invariant mass of the K~^tt~ fj.~ system vs. (f' for signal candidates after 
the full selection, where the )• signal decay is clearly visible as vertical band. The 

regions 8.0 < < ll.OGeV^/c'^ and 12.5 < < 15.0GeV^/c^ contain the tree-level decays 

B^ ^ J/V’(—)• and V'(25)(—which are used as important control 

decays but vetoed when selecting B^ —)• signal candidates. Integrated over q'^, the 

B^ ^ signal yield is 2 398 ± 57, as shown in Fig. 1. 

2.3 Angular analysis 

The angular observables are determined by performing an unbinned maximum likelihood fit in 
bins of q‘^, using a q^ binning which is both finer and more regular than the binning used in 
Refs. The ht uses the invariant mass distribution of the K'^tt~ ij,~ system, the invariant 
mass and the three-decay angles as input distributions without applying angular foldings. 
This allows to quote the covariance matrices for all eight angular observables which is important 
for the use of the results in global fits. The invariant iF+7r“ mass distribution is used to constrain 
the contribution from events where the K^'k~ system is in a spin-0 configuration, the so-called 
S-wave. The additional six parameters for the description of the S-wave and the interference 
terms with the K*^ P-wave are treated as nuisance parameters and allowed to vary in the ht. The 
trigger, reconstruction and selection of signal events causes distortions of the distributions of q^ 
and the decay angles. This acceptance effect is modelled using a multidimensional combination 
of Legendre polynomials. The polynomial coefficients are determined using a moments analysis 

































Figure 2 - The angular observables Fl, Afb, Ss and P 5 overlaid with SM predictions from (purple) Ref. and 
(orange) Ref. . 


of a large sample of simulated —)• events, generated according to a phase-space 

model. The Feldman-Cousins method is used to guarantee correct coverage for the angular 
observables even for low signal yields. 


2.4 Results 

The results for Tp, ^fBj S 5 and P 5 are given in Fig. 2, overlaid with SM predictions from Refs.^’^^ 
and Ref. The longitudinal polarization fraction Fp and the forward-backward asymmetry 
Afb are found to be in good agreement with SM predictions. Interestingly, for ^fBj the data 
points seem to lie systematically below the predictions in the 1.1 < < 6 GeV^/c^ range. The 

measurement of the less form-factor dependent observable P 5 is found to be compatible with 
the previous publication”^ and lies above the SM prediction in the 4.0 < < 8.0GeV^/c^ 

region. The deviation from the SM prediction corresponds to 2.9 cr for each of the two q^ bins in 
the region 4.0 < q^ < 6.0GeV^/c^ and 6.0 < q^ < 8.0GeV^/c^. Neglecting correlations between 
the bins, the probability to find a deviation of this size or larger for two degrees of freedom 
results in a naive significance of 3.7cr. The remaining observables S 3 , S 4 , Sj, Sg and Sg are 
given in Ref. and show good agreement with SM predictions. 


3 Branching fraction measnrements of R—> and decays 

Gompared to angular observables, branching fraction measurements of 6 —?• processes 

tend to have larger associated theory uncertainties, since they are directly impacted by the 
hadronic form-factors. However, the impact of theory uncertainties can be mitigated, by per¬ 
forming measurements of ratios of branching fractions where form-factor uncertainties cancel at 
leading order. Examples of such quantities are the isospin asymmetry Aj and the CP-asymmetry 
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Figure 3 - Differential branching fraction of (left) —>• fi (middle) B° —>■ and (right) B'^ —> 

K*'^, overlaid with SM predictions 
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In Refs.^®’^®, Aj and ^cp are found to be compatible with SM predictions^^’®. The corresponding 
differential branching fraction measurements for the rare decays R"*"—)• fi'^ii~, B^^ 

and R'*'—>■ K *~^are given in Fig. 3. They are compatible with, but tend to lie below, SM 
predictions 

Using lfb~^ of data taken during 2011, LHCb also determines the differential branching 
fractions for the rare decays )• and B^^ The differential branching 

fractions tend to be below SM predictions both at low q^, where updated light cone sum rule 
calculations are available and at high (f‘, where lattice calculations exist For the decay 

B^ —>■ the tension in the region 1 < < 6 GeV^/c^ corresponds to 3.1 a. It is interesting 

to note, that the deviation of the branching fractions points to a deviation of the b —)• 
couplings which is compatible with, but less significant than, what is observed from the angular 
observables in )• at low Updated measurements of B{B^^ ^~) and 

—)• using the full Run I data sample are currently in preparation to clarify the 

situation. 


4 Branching fraction of B^^^ tt+tt 

The TT~^TT ~final state can be reached from both the decay of a B^ meson and the decay 
of a Bg meson. The B^ decay is expected to be dominated by the b —)• transition 

BO ^ the decay by the b —)• spApT transition —)• p~. While b ^ d 

decays are expected to be suppressed by the factor |VtdP/|ktsP ~ 0.04 compared to 6 —)• s 
transitions in the SM, this is not necessarily the case for SM extensions. 

The TT~^TT~p'^p~ final state is studied using the full Run I data sample taken by the LHCb 
experiment®^, corresponding to an integrated luminosity of 3fb~^. The invariant mass of the 
7 r+ 7 r“ system is required to be in the range 0.5 — 1.3GeV/c^ containing both the as well as 
the f^{980) resonance. Figure 4 gives the invariant mass distribution of the p~^p~ system 
for the charmonium modes —)■ J/'i/) 7 r^ 7 r“, that are used as control decays for the fit model, 
as well as the signal decays B^^-^ —>■ 'k~ p~^p~. The signal yields are found to be 40 ± 10 ± 3 
for the B^ —)• p~ decay and 55 ± 10 ± 5 for the —)• TT~^7r~p'^p~ decay, resulting 

in significances of 4 . 8(7 and 7.2 a, respectively. The branching fractions are determined with 
respect to the normalisation mode B^ ^ J/il^K*^. They are found to be 

^ 7r+7r-//+;U-) = (8.6 ± 1.5,tat. ± 0.7syst. ± O.Tnorm.) X 10-®, 

^ TT+TT-p+p-) = (2.11 ± 0.51stat. ± 0.15syst. ± 0.16norm.) X 10“®, 









Figure 4 - Invariant mass of the tv'^'k /i final state for (left) tree-level charmonium decays 
and -Bs J/'ip'K'^'K~ and (right) the rare decays >■ 7r"'"7r“/r"'"^“ and -8°—>■ iv'^ n~ fi~. 
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Figure 5 - The (left) differential branching fraction and (middle) leptonic, as well as the (right) hadronic forward- 
backward asymmetry, overlaid with SM predictions 


when correcting for the (f' regions removed by the vetoes of the and 

decays, in agreement with SM predictions 

5 The rare baryonic decay —)• 

The study of the rare decay A^ —)• is of particular interest due to the half-integer 

spin of the A^ baryon and the hadronic dynamics involving the heavy b and a light diquark 
system. Furthermore, the A decays weakly into the p'K~ final state, allowing access to new and 
complementary information compared to mesonic b —?■ decays^®. 

The decay was previously studied in where no evidence for signal in the region 
below the J/'0 was found. An updated analysis is performed, using the full LHCb Run I data 
sample Figure 5 gives the differential branching fraction. Evidence for signal is found below 
the charmonium resonances at low q^, the differential branching fraction for the high range 
15 < < 20GeV^/c'^ is determined to be (1.18^Q;g8 ± 0.03 ± 0.27) x 10“'^ GeV^c^^^. Angular 

analyses are performed for the q'^ bins where evidence for signal is found and the angular ob¬ 
servables Apg and Apg, the forward-backward asymmetries in the dimuon and hadron system, 
are determined. As shown in Fig. 5, Apg and Apg are found to be in good agreement with SM 
predictions ^^, 44 ^ 


6 A test of lepton universality using the decay K+e+e 

The ratio Rk in the region g'^ax] is defined as 

ldg2 


Rk = 


flmax 
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^ . 
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dT[B+^K+e+e- 


-dq^ 
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where F denotes the q^ dependent partial width. Due to the universal coupling of 7 and Zq 
bosons to leptons, Rk in the region 1 < < 6 GeV^/c^ is predicted to be one with an uncertainty 

















































xlO’ 




-»-LHCb -B-BaBar -*-Bclle 



LHCb : 


1 


1 , 

t 

SM 




O' . 

0 5 10 15 20 


q- [GeV^/e^] 


Figure 6 - (Left) —>■ J/tpK'^ and (middle) S'*" —>■ K^e^e signal candidates, triggered on the electron in the 
final state. (Right) The ratio Rk as determined by LHCb"^^, BaBar^* and Belle for different ranges. 


of less than 10 ^ 45,46^ Small corrections to the ratio arise only from phase-space effects and 
Higgs penguin contributions. 

The measurement is experimentally challenging due to a lower trigger efficiency for electrons 
compared to muons and the higher emission of Bremsstrahlung which deteriorates the resolution 
of the invariant mass of the K'^e^e~ system. In the range 1 < < 6 GeV^/c^, Rk is determined 

to be 


Rk = 0.745^q Q 74 (stat.) ± 0.036(syst.), 

which corresponds to a deviation of 2.6 u from the SM predictionFigure 6 shows the decay 
B'^ —)■ J/'i/i K~^ , which is used to study the effect of Bremsstrahlung and understand the relative ef- 
hciency between reconstructing dimuon and dielectron modes. The signal decay B'*'—)• K~^e'^e~ 
is also given, as well as the LHCb measurement of Rk in comparison with results from the 
B factories Since Rk is free from hadronic uncertainties, the result received considerable 
attention from theory Further tests of lepton universality are in preparation, 

including the measurements of Rk* and R^. 

7 Angular analysis of ^ K*^e'^e~ 

The study of rare decays with electrons in the final state allows to perform analyses at very 
low q^, due to the tiny electron mass. At low g^, the contribution from Feynman diagrams in 
which a virtual photon couples to the lepton pair dominates. This allows to probe the photon 
polarisation, which is left-handed in the SM. 

LHCb performs an angular analysis of the decay B^ —>• K*^e'^e~ in the g^ range 0.002 < 
g^ < 1.120 GeV^/c^ The four angular observables Fl, A^\ and are determined 
from an unbinned maximum likelihood fit to the decay angles cos 9i, cos 6k and (f>. Of particular 
interest are the observables A^^ and A^ that are sensitive to the photon polarization. Figure 7 
gives the angular fit projections. The measured angular observables are 

Fl = -b0.16± 0.06 ±0.03 
A^^^ = -0.23 ± 0.23 ± 0.05 
= ±0.10 ±0.18 ±0.05 
A^ = ±0.14 ±0.22 ±0.05, 

which is in good agreement with SM predictions®^’^®. The constraints from AI^ on the 
contributions from right-handed currents are more precise than those obtained from the average 
of the time dependent CP-asymmetries in radiative )• ?• decays®®’®®. 




















Figure 7 - The three decay angles (left) cos 9i (middle) cos6k and (right) (p of >■ K*°e^e signal candidates, 
overlaid with the projections of the htted probability density function. 


8 Conclusions 

Most of the observables in rare decays are found to be in good agreement with SM predictions. 
However, three interesting tensions emerge: An update of the angular analysis of the decay 
BO j^*o 

confirms a deviation of the angular observable in the two bins 4 < 

< 6GeV^/c^ and 6 < < 8GeV^/c^, with a significance of 2.9 ct in each; Furthermore, the 

branching fraction of the rare decay Bg —)• in the range 1 < < 6 GeV^/c^ is 3.1 cr lower 

than a recently updated theory prediction; Finally, the measurement of Rk shows a tension 
with lepton universality at 2.6 cr. 

Consistent NP explanations of all observed tensions in rare decays exist, and first global fits 
including the updated results on —)• angular observables prefer the NP solution 

over the SM by 3.7 . However, it is too early to speak of clear signs of processes beyond the 

SM; Unexpectedly large hadronic contributions still can not be excluded. The results clearly 
motivate future work, both from theory, as well as from experiment, where complementary 
measurements of rare b —)• (s, d)i£ processes will be performed. For the Run I LHCb data, this 
includes an update of the analysis of the decay Bg —)• (j)p ^and an updated branching fraction 
measurement of the decay )• K*^fj.~. In addition, further tests of lepton universality and 
lepton number violation are in preparation. The data sample LHCb will collect during Run H 
will further improve the experimental sensitivity and allow to probe the operator structure of 
rare decays with unprecedented precision. 
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